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A TEST OF THE CENTRAL-MARGINAL MODEL 
USING SAND SCORPION POPULATIONS 
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Tsunemi Yamashita 1 and Gary A. Polis: Department of Biology, Box 1812 Station B, 
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ABSTRACT. The central-marginal model proposes marginal populations contain a lower density of individ¬ 
uals, lower levels of genetic variation, and are more isolated than populations in the center of a species range. 
Previous tests of the model used Drosophila , organisms capable of extended dispersal. We test the central- 
marginal model with scorpions, organisms with restricted dispersal abilities. We measured genetic variation 
through allozyme analysis of eight loci (five polymorphic, three monomorphic) to obtain estimates of hetero¬ 
zygosity. We compared differences between the two types of populations with a split-plot ANOVA. We also 
compared central and marginal populations using standard parametric tests. We found marginal populations 
contain lower genetic variation than central populations. These populations may be important as models in 
conservation to study the effects of fragmentation. 


Peripheral or marginal populations are those 
on the boundaries of a species’ geographic range. 
They exhibit unique properties not evident in 
populations in the center (Brussard 1984). Gen¬ 
erally, as one moves outward from the center of 
a species’ range, populations are hypothesized to 
become less dense, more isolated, and less vari¬ 
able genetically within populations (da Cunha et 
al. 1959; da Cunha & Dobzhansky 1954; Carson 
1959; Soule 1973; Brussard 1984). These trends 
are embodied in the central-marginal model (Le- 
wontin 1974). Several explanations exist. One 
hypothesis (da Cunha & Dobzhansky 1954) states 
genetic polymorphism is positively correlated 
with the number of niches an organism occupies 
— more niches are available at the center of a 
species range. Yet Brussard’s (1984) free recom¬ 
bination hypothesis proposes linkage disequilib¬ 
rium is favored in central populations because 
extreme phenotypes are selected against, i. e., 
stabilizing selection occurs. Here, favored genes 
are linked together to create a stable phenotype 
with maximal fitness. In theory, at the range mar¬ 
gins, linkage equilibrium is selected in stressful 
environments to create novel phenotypes better 
able to survive sub-optimal conditions (Brussard 
1984). 

Research with Drosophila suggests the predic- 
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tion that allozyme heterozygosity is reduced at 
range margins is not valid (Brussard 1984). How¬ 
ever, Drosophila probably disperse relatively great 
distances compared to most non-flying taxa. Non¬ 
flying taxa (vertebrates) show a decline in allelic 
diversity in marginal populations (Soule’ 1973). 
We propose scorpions, unable to disperse great 
distances, also may illustrate the predictions of 
Lewontin’s (1974) central-marginal model. 

METHODS 

Relevant scorpion biology.—The scorpion Pa- 
ruroctonus mesaensis Stahnke is restricted to sand 
dunes and sandy substrates scattered throughout 
the southwestern United States and northern 
Mexico. It occurs in dense populations (range = 
1600-5000/ha) (Polis & Yamashita 1991) and it 
is an ecologically important species, as a gener¬ 
alist predator in desert food webs (Polis 1979; 
Polis & McCormick 1986, 1987). This scorpion 
is dispersal limited because specialized morpho¬ 
logical features adapt and restrict it to sand. The 
species possess numerous modified setae on their 
tarsi (sand shoes) to facilitate sand movement 
and burrow construction (Polis et al. 1986). Fur¬ 
ther, their ability to detect substrate vibrations 
to localize prey only functions well on sand 
(Brownell & Farley 1979). These specializations 
reduce the likelihood of extended migration. 
However, reproductive males move extensively 
when searching for mates (Polis & Farley 1979, 
1980). 
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Collection sites: 

Central Populations: 

A. Glamis 

B. Yuma Desert 

C. Yuma Army 

D. Ogilby 

E. N Fontana 

F. Glamis East 


Marginal Populations: 

G. Windy Point 

H. Kelso 

I. Needles 

J. Bouse 

K. Phoenix 

L. Agua Caliente 


0 40 
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Figure 1 .-A partial map of the desert areas of California, Arizona, and Baja California Norte where populations 
of Paruroctonus mesaensis were collected. A proposed range of P. mesaensis is indicated by the dashed line. 


Collection and electrophoresis procedures.— 
We collected the scorpions over three years, 1989- 
1991. Paruroctonus mesaensis is easy to collect 
as it (and all scorpions) fluoresces when illumi¬ 
nated with a ultraviolet light (Sylvania F8T5/ 
BLB). We transported the scorpions alive to the 
laboratory, froze them with liquid nitrogen, and 
stored them at —70 °C until electrophoretic anal¬ 
ysis. 

We electrophoresed five polymorphic and three 
monomorphic loci from 28 populations. This re¬ 
search was the first allozyme analysis of any scor¬ 
pion populations, and these eight loci were the 
only resolvable loci from a screening of 25 en¬ 
zyme loci on 11 different buffer systems. We as¬ 
sayed a mean of 29 individuals (± 5.7) from each 
population for each locus. The specific protocols 


and other pertinent methodology are described 
in Yamashita (1993). 

Central and marginal population determina¬ 
tion.— We determined the geographical center 
of P. mesaensis populations from collection data 
(Haradon 1983; D. Gaffin pers. comm.; Yama¬ 
shita 1993) (Fig. 1). We determined the following 
range extremes: in the north (Death Valley, Cal¬ 
ifornia); in the south, Cabo Lobos (Sonora, Mex¬ 
ico); in the east (Phoenix, Arizona); and in the 
west (Windy Point in the Coachella Valley, Cal¬ 
ifornia). We determined the range center to be 
located 65 km north of Yuma, Arizona by lo¬ 
cating the midpoint between the extreme north¬ 
ern and southern and the western and eastern 
populations. We designated six populations 
within a 32 km radius of the range center as 
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Central populations: Glamis, Yuma desert, Yuma 
Army, Ogilby, N. Fontana, and Glamis East. We 
delineated peripheral populations as those near¬ 
est to edges of the range. These include Windy 
Point, Kelso, Needles, Bouse, Phoenix, and Agua 
Caliente. 

Statistical analyses.—We performed two types 
of analyses to determine if differences between 
central and marginal populations exist. The first, 
a method outlined by Weir (1990), tested if het¬ 
erozygosity differences exist between population 
types. This design, similar to a split-plot ANO- 
VA, considered variation from five sources: pop¬ 
ulations, individuals within populations, loci, loci 
by populations, and loci by individuals within 
populations. Heterozygotes are entered as l’s, 
homozygotes, as 0’s and the data from each cen¬ 
tral or marginal populations were pooled. We 
used electrophoretic data from 20 individuals 
from each of six central and six marginal pop¬ 
ulations (120 individuals in each population type) 
and five loci. We used Weir’s analysis because 
estimates of heterozygosity often exhibit large 
interlocus variances and non-normal distribu¬ 
tions; therefore, many standard parametric tests 
may be inappropriate (Archie 1985). 

Second, we performed standard parametric 
statistics to determine differences between cen¬ 
tral and marginal populations. The variables ex¬ 
amined were observed average heterozygosity, 
mean allele number, and percent polymorphism 
of each population. We calculated these variables 
using BIOSYS-1 (Swofford & Selander 1989). 
Observed heterozygosity per locus is the fraction 
of heterozygous individuals from a given sample 
for a particular locus (Ferguson 1980; Weir 1990). 
Observed average heterozygosity is the mean 
value from all loci. Although heterozygosity val¬ 
ues commonly undergo an arcsine transforma¬ 
tion, our data did not require such a procedure 
because most heterozygosity values fell between 
0.30 and 0.70, a range that does not require trans¬ 
formation (Sokal & Rohlf 1981). 

We used two other indices of genetic vari¬ 
ability. Percent polymorphism is the mean num¬ 
ber of polymorphic loci in a population. Here, a 
locus is polymorphic if the frequency of the most 
common allele is 0.95 or less. The mean number 
of alleles per locus is the number of alleles at 
each locus averaged across all loci. 

RESULTS 

The split-plot ANOVA and standard para¬ 
metric tests showed the mean heterozygosities 


Table L—Split plot ANOVA analysis of two pop¬ 
ulation categories (Central and Marginal). This design 
is taken after Weir 1990. 



df 

MS 

F-value 

P-value 

Category 

Individuals within 

1 

3.54 

21.60 

<0.001 

categories 

222 

0.166 

1.01 

>0.437 

Loci 

4 

4.29 

26.17 

<0.001 

Category x loci 

4 

1.29 

7.86 

<0.001 

Error 

888 

0.164 




from the central and marginal populations were 
significantly different (Tables 1, 2). The results 
of the split-plot ANOVA (Table 1) show values 
from all levels of analysis (populations; individ¬ 
uals within populations; loci; population x loci) 
were significant (P < 0.001) except individuals 
within populations (P > 0.437). These results 
establish that central populations are signifi¬ 
cantly different from marginal ones in hetero¬ 
zygosity. Furthermore, the significant among loci 
effects suggests that each locus expressed a dif¬ 
ferent pattern of heterozygosity from other loci. 
Loci within marginal populations were signifi¬ 
cantly different from loci within central popu¬ 
lations, which suggests within each population 
type (central or marginal), the same locus ex¬ 
pressed significantly different heterozygosities. 

The mean genetic variability (observed het¬ 
erozygosity) in marginal populations (0.106 ± 
0.025, n = 6) was significantly less than central 
populations (0.164 ± 0.018, n = 6; t = 3.80,0.05 
> P > 0.01) supporting the central-marginal 
model. The mean allele number for the central 
populations (1.73 ± 0.144) was marginally great¬ 
er (t = 2.21, 0.1 > P > 0.05) than that of mar¬ 
ginal populations (1.52 ± 0.095). The mean per¬ 
cent polymorphism for central populations (45.83 
± 6.45) was also marginally greater (35.42 ± 
9.41, t = 2.71,0.05 > P > 0.01). 

DISCUSSION 

The significant differences between central and 
marginal populations for allozyme heterozygos¬ 
ity, mean allele number, and percent polymor¬ 
phism are consistent with the central-marginal 
model (Brussard 1984). Paruroctonus mesaensis 
is one of the species that fits the predictions of 
this model; tests of the model using Drosophila 
allozymes failed to exhibit similar patterns 
(Brussard 1984). In our study, the decrease in 
genetic variability in marginal populations prob- 
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Table 2.—A comparison of genetic variability between central and marginal populations. See text for discus¬ 
sion. 



Mean 

heterozygosity 

Mean allele 
number 

% Polymorphism 

Central populations 

Glamis 

0.153 

1.75 

37.5 

Yuma Desert 

0.176 

2.00 

50.0 

Yuma Army 

0.136 

1.63 

50.0 

Ogilby 

0.167 

1.63 

37.5 

N. Fontana 

0.166 

1.63 

50.0 

Glamis East 

0.187 

1.75 

50.0 

Mean 

0.164 

1.73 

45.8 

SD 

0.018 

0.144 

6.5 

Marginal populations 

Windy Point 

0.103 

1.50 

25.0 

Phoenix 

0.078 

1.38 

25.0 

Bouse 

0.131 

1.63 

50.0 

Agua Caliente 

0.140 

1.50 

37.5 

Kelso 

0.084 

1.63 

37.5 

Needles 

0.097 

1.50 

37.5 

Mean 

0.106 

1.52 

35.4 

SD 

0.025 

0.095 

9.4 

Central vs marginal 

/-statistic 

3.80 

2.21 

2.71 

P values 

0.05 > P > 0.01 

0.1 > P > 0.05 

0.05 > P > 0.01 


ably stems from reduced gene flow or smaller 
overall population size. Because scorpion dis¬ 
persal is local, populations at the range margin 
are less likely to receive migrants from other pop¬ 
ulations compared to more central populations. 
Central populations exhibit the highest allele 
number and percent polymorphism. These pop¬ 
ulations may maintain higher genetic variability 
because exchange with other nearby populations 
is more frequent and population size is generally 
larger in the center of the range. However, mod¬ 
els suggest that a very small effective population 
size (n e < 10 individuals) is required to reduce 
significantly the number of alleles per locus with¬ 
in a population (Nei et al. 1975; Rice & Mack 
1991). 

Some marginal populations (Needles, Bouse, 
and Phoenix) are geographically isolated from 
other populations. Needles, north of a mountain 
range present on either side of the Colorado Riv¬ 
er, is effectively isolated. The Bouse population 
exists on the eastern edge of the Cactus Plain, a 
large sandy region in western Arizona. It is sur¬ 
rounded by rocky habitat and isolated from the 
nearest population by 40 km. Although Bouse is 
not separated by a large distance, the interme¬ 
diate rocky substrate effectively curtails dispers¬ 


al. No P. mesaensis were observed in > 30 hours 
of searching on rocky habitats adjacent to sandy 
areas (Polis, unpubl. data). 

The Phoenix population is the most eastern 
and one of the most genetically depauperate pop¬ 
ulations. Emigration into this area probably oc¬ 
curred along the dry river beds of the Salt and 
Gila rivers. This population relies on unidirec¬ 
tional gene flow since the substrate outside the 
river bed is a dispersal barrier to the psammo- 
philic scorpion and no populations exist to the 
east. The low values of genetic variability may 
be a result of two primary factors: it is a periph¬ 
eral population with a low population size and 
receives little gene flow from other populations. 

Peripheral isolate formation may have been 
enhanced by large scale floods in the Holocene 
(Ely et al. 1993). The Salt and Gila rivers ex¬ 
perienced large-scale floods in the last 5600 years 
(Ely et al. 1993). Periods of minor and major 
floods were interspersed. These extreme floods 
may have created expansion corridors for scor¬ 
pion movement and isolated populations by 
fragmenting previous habitat. 

Although several studies have compared cen¬ 
tral and marginal populations, a clear trend is 
not evident (Hoffman & Parsons 1991). Some 
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report no decrease in genetic variability within 
marginal populations, e, g., Drosophila alio- 
zymes (Bmssard 1984) and an annual grass (Rice 
& Mack 1991), Analysis of dispersal limited an¬ 
imals (frogs) report a decrease in genetic vari¬ 
ability among marginal populations (Sjogren 
1991). Further research into the properties of 
marginal populations is warranted because mar¬ 
ginal populations, with their insular or penin¬ 
sular properties, are similar to populations frag¬ 
mented through man's encroachment upon the 
environment. 
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